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Abstract 

The presence of a bulk viscosity for the cosmic fluid on a single Randall-Sundrum brane is 
considered. The spatial curvature is assumed to be zero. The five-dimensional Friedmann equation 
is derived, together with the energy conservation equation for the viscous fluid. These governing 
equations are solved for some special cases: (i) in the low-energy limit when the matter energy 
density is small compared with brane tension; (ii) for a matter-dominated universe, and (iii) for a 
radiation-dominated universe. Rough numerical estimates, for the extreme case when the universe 
is at its Planck time, indicate that the viscous effect can be signiflcant. 
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I. INTRODUCTION 

As is known, the proposal of Randall and Sundrum (RS) Q to solve the hierarchy problem 
- in contrast to the approach of Arkani-Hamed et al. jsj, for instance - was to introduce 
a warped metric in five-dimensional space, the fifth dimension y being compactified on an 
orbifold /Z2 of radius R such that —-kR < y < ttR. In the two-brane model (RSI) the 
fixed points y = and y = ttR were the locations of the two three-branes, whereas in the 
one-brane model (RS2) the brane was taken to be located at y = 0. 

In a more realistic physical setting where the cosmic fluid is assumed to be present on 
the brane, the original RS warped metric has to be generalized. In addition to the brane 
tension (or four-dimensional cosmological constant), which we will call a, one has to allow 
for a fluid pressure p and a mass-energy density p in physical space. The form of metric 
that we shall consider in the following to describe this situation, is 

ds"^ = -n^{t, y)dt^ + a^{t, y)6ijdx'dx^ + dy"^, (1) 

where n(t, y) and a(t, y) are functions to be determined from Einstein's equations. This is 
a kind of metric that has been made use of earlier - cf. Refs. [3-6], for instance - but we 
shall here specialize to the case of zero spatial curvature, /c = 0. There is by now ample 
experimental evidence supporting the flat space assumption. 

We will in the following analyze a single flat membrane situated at the position y = 
in such a flat space. An isotropic viscous fluid with bulk viscosity C is taken to be present 
on the brane. Now, there are in general two different viscosity coefficients for a fluid, but 
we will henceforth ignore the second component, i.e., the shear viscosity rj. The reason for 
this is the assumed perfect isotropy of the fluid. Actually this is a rather delicate point, 
due to the much larger magnitude of rj than of C,. Thus if one considers the conventional 
plasma epoch in the history of the universe a few tenths of seconds after the big bang, when 
the viscosity coefficients are calculable from ordinary kinetic theory, then one finds that 77 is 
larger than C, by more than 15 orders of magnitude. Even a small amount of anisotropy in 
the fluid would thus easily outweigh the effect of the bulk viscosity. This point is discussed 
more closely in Refs. 0] and Q]. 

The case of a bulk-viscous cosmic fluid has to our knowledge not been much studied 
studied before in a five-dimensional context, although we have recently become aware of 
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the recent papers of Harko et al. in this direction j9|, llO|. Some of our own motivations for 
undertaking this kind of study were the following: 

First, there is the obvious need of generalization of the theory. The case of a nonviscous 
fluid is from a hydrodynamical point of view a very idealized situation, although it is a 
useful approximation in many cases. A condition for five- dimensional viscous cosmology to 
have physical meaning, is that it is robust enough to encompass the presence of viscosity 
coefficients. Introduction of these coefficients implies that we are accounting for thermal 
irreversibility to the first order in the fluid velocity's derivatives. 

Our second point relates to entropy. The nondimensional entropy per baryon in the 
universe is known to be very large, of the order of 10^ ^|. Although this entropy is 
basically to be associated with microscopical processes one would expect, by analogy with 
ordinary hydrodynamics, that it is describable in terms of phenomenological parameters 
such as viscosity. 

Finally, there are the energy conditions, the most important of which is probably the 
weak energy condition (WEC), p + p >0 The extent to which the presence of viscosity 
influences the energy conditions seems to be to a large extent unexplored. Some discussion 
in the case of anisotropic cosmology can be found in [13]. The energy conditions are likely 
to be of importance also for the recent theories about phantom matter [l^ . 



II. BASIC CONCEPTS 

We summarize our assumptions: The metric, as mentioned, is given by Eq. (^, describing 
the space on the brane itself {y = 0) as well as in the surrounding bulk, the latter being 
empty except from the presence of a flve- dimensional cosmological constant A. On the brane, 
the tension a is also assumed to be constant. The metric extends out to the horizon, 
determined by goo = 0. The flve-dimensional Einstein equations are 

Rmn — -9mnR + 9m = /^"^Tmn, (2) 

where capital Latin indices denote flve-dimensional indices {t,x'^,x'^,x^,y), and = 871G5 
is the flve-dimensional gravitational coupling. The energy-momentum tensor describes every 
classical source except from the "vacuum" fluid, i.e., A. 

The components of the Ricci tensor are conveniently found using the Cartan formalism. 
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Since the relationship between the components T^.jfj of the energy-momentum tensor in 
orthonormal basis and the components Tmn in coordinate basis are 



J-tt = n Iff, lij = a ifj, J-ty = ^J-iyy ''■yy ~ ^yy^ 



(7) 



we get from Eq. Q 




a" a'^ 
a \ a . 



An^ = K^Ttu 



(8) 



5ij < a^ 



a' fa' ^ 2n'" 
a \ a n 



\ 2a" n". 


a" 


d 


( a 


2h\ 


2dn 


+ — + — 


+ — 






+ — - 




/ a n J 




■a 


\ a 


n J 


a ■ 



g / d ra' b! 
\a n a , 



tyi 



3a' ( a' n' 
— - + - 

a \ a n 



+ 



A = /t^T 



yy 



«'T,,, (9) 



(10) 
(11) 



a /a n 
n^ la \a nj a\ 

Here a = a{t, y), n = n{t, y), overdots and primes meaning derivatives with respect to t and 
y repectively. So far, the form of the energy-momentum tensor has not been specified. 

A note on dimensions: As the spatial curvature A; = 0, it is convenient to let the di- 
mensions be carried by the coordinates so that a and n become nondimensional. Thus 

[a] = [n] = 1, [x^] = [y] = cm, [A] = cm^^, [k^] = cm^. 



III. ENERGY-MOMENTUM TENSOR AND GOVERNING EQUATIONS 

The governing equations for the metric components are obtained by integrating Einstein's 
equations across the singular boundary at ?/ = 0. As for the energy-momentum tensor, 
the contribution to it from the constant brane tension is Tmn = ~^{y)<^ ^m^n d^f^- 
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orthonormal frame this means that T^^ = S{y)a, T-f- = —S{y)cr, which effectively imphes 
j> = —a on the brane, i.e., the conventional pressure-energy relationship for a domain wall 

Adding the energy-momentum tensor T^""^ for the viscous fluid, we get 



Tmn = 5{y){-ag,, + T^'^Wn. (12) 



where we insert the conventional four-dimensional expression 

TX'' = pU,U, + {p-CO)K,. (13) 

Here = (f/°, W) is the fluid's four-velocity, 6 = U^-^^ is the scalar expansion, and 
h^v = Qfiu + U^Uy is the projection tensor. We work henceforth in an orthonormal frame, cor- 
responding to = (1, 0, 0, 0). With the notation ao(t) = a(t, y = 0) and no{t) = n{t, y = 0) 
we can write the metric on the brane as 

ds"^ = -nl{t)dt^ + al{t)6ijdx'dx\ (14) 

The scalar expansion is calculated as ^ = 3do/ao + ho/nQ. We impose the gauge condition 
no(t) = 1 which physically means that the proper time on the brane is taken as the time 
coordinate, and so get 

e=^-^^ 3Ho. (15) 

This is the same relationship as in ordinary FRW cosmology. Letting p denote the effective 
pressure, 

p = p- 3HoC, (16) 

we thus have 

Tr" = p, 7;f^=pa^5,,. (17) 

Consider next the junction conditions. The metric is continuous across the brane, but its 
derivatives are not. The distributional parts of a" and n" have to be matched with the 
distributional parts of the energy-momentum tensor. Following Ref. P| we write a" = 
d" + [a']5(y), where [a'] = a\y = 0"*") — a'{y = 0~) is the jump across y = and d" is the 
nondistributional part. Similarly, we write n" = n" + [n']6{y) . The matching procedure 
applied to Eqs. (jH} and yields 

^^ = -l^\a + pl (18) 
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for the nondistributional parts. (For simphcity we have written a" and n" instead of a" and 
h" in Eqs. ()20|) and ((21]) •) In the junction conditions (fTHj) and (fT^ we have made use of 
hq = 1. Equation ()18|) is seen to be formally the same as for an ideal fluid, whereas Eq. ()19|) 
is explicitly viscosity dependent. 

We assume henceforth that there is no flux of energy in the y direction. It means. 



Tty = 0. 
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Equation (fTUI) thus yields 
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for arbitrary y. 
By means of Eq. 



we can write Eq. (p])) as 



d 

dy 



K)^-(v)' 



2 d 1 4 

— A— -a 

3 dy% ^ 



Integrating this equation with respect to y we get 
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where C = C{t) is an integration constant. 

Assuming Z2 symmetry y — > —y for the scale factor a across the brane as well as continuity 
for y —>■ 0, we obtain from the boundary condition (|TH|) 

1 



(27) 



^J = A + -.Vp+-.V + ^ (28) 



at y = 0. Also by evaluating Eq. (26) on the brane and removing the y derivatives by 
inserting Eq. ()27|) we finally obtain (remembering no{t) = 1) 

(recall that subscript zero refers to the brane). We have here introduced the quantity 

A = -A+— «;V, (29) 
6 36 ' ^ ^ 

which can be interpreted as the effective four-dimensional cosmological constant in the five- 
dimensional theory. The "new" Friedmann equation (j28|) can be contrasted with the con- 
ventional Friedmann equation in four-dimensional cosmology (when k = 0): 

Ho = l^A + l4p^ (30) 

where A4 and k| = 87rG4 are four-dimensional quantities. The essential new properties of 
Eq. (PHj) are thus the presence of a p^-term, and also the presence of a "radiation" term 
C/al. 

The (nondistributional) Einstein equation (j^UI) can be solved analytically for a{t,y). In- 
serting Eq. into Eq. (pUj) and multiplying with we obtain 

{aa'Y + -Aa^ -dl = 0. (31) 
3 

Let us first assume that A > 0, i.e., the de Sitter (dS) case, and let us take as trial solution 

a^{t,y) = Acos{2fidy) + B sm{2i^dy) + D, (32) 



where A, B, D are arbitrary functions of time, and where /x^ = y A/6. In order for Eq. 
to satisfy Eq. ()31|). we must take 
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where the elimination of do in the final step is done by means of Eq. ()28j) . 

The functions A and B in Eq. (p?^ are determined by using the assumption of Z2 symme- 
try, together with the junction condition (fTHj) . Writing = A\y=Q+, A~ = A\y=Q-, B^ = 
-B|y=o+ and B^ = i?|j^=o-) ^^id using simple symmetry relations for the trigonometric func- 
tions, we obtain from Eq. ()32j) 

A+ = A~ = A, 5+ = -B- = B. (34) 



After a simple calculation using Eqs. (fTHj). (jH^ . and (jH^ . we find 
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The full expression for a(t, y) is thus 
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The corresponding solution of Einstein's equation for an AdS bulk space (i.e., A < 0) is 
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(38) 



where /x = y — A/6. Equations (fHTjl and (jHEl) are the same solutions as given in Ref. j^, but 
here with a fiuid on the brane (i.e., cr — > a + p). 

The second junction condition, Eq. ()19p. relates to the conservation equation for energy. 
Using the same trial solution as before, i.e., Eq. (jH^ for a dS bulk space (which we know 
satifies the Einstein equation as long as D is given by Eq. ()33p). and taking into account 
the relation 

d\t,y) 
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which follows from Eq. by different at ion, we find by evaluating n' on both sides of the 
brane 

p + 3{p + p)Ho-9CH^,=0. (40) 

The calculation is analogous to that leading to Eq. but now with the junction condition 
for n, Eq. p9j) . taken into account. 



8 



■ liquation ()4Up is the same as that obtained in ordinary four-dimensional cosmology 



16l . llTj l . This is somewhat surprising, as there is seemingly no simple physical reason why it 
should be so. 

The anti-de Sitter (AdS) universe, A < 0, yields the same result. 



IV. COSMOLOGICAL SOLUTIONS 



We introduce the equation of state in the conventional form 

P = (7 - (41) 

where 7 is a constant lying in the interval < 7 < 2. As from now on we will consider 
cosmological solutions only, we will henceforth leave out subscripts zero referring to the 
brane position y = 0. 

As is commonly believed, the evolution of our universe consists of four epochs: (1) the 
early inflationary epoch, (2) the radiation dominated epoch, (3) the matter dominated epoch, 
and (4) the present mini- 

energy density of our universe is dominated by the effective cosmological constant in epochs 
(1) and (4). We will assume that the influence of the radiation term C/a^ is negligible except 
in the radiation era, epoch (2). 

With C = the governing equations reduce to the set 

H' = \ + ^K'ap+^K'p', (42) 

p + ^fpO-ce^ = 0, (43) 

with e = 3H. 

Before embarking on the viscous case, let us briefly consider the situation when C = 0. 



A. Nonviscous fluid 



From Eqs. ()42j) and ()43j] we get for the cosmological time 

dx 



(44) 
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where x = p/o" is a nondimensional parameter. From Eq. ()43p it follows that when C = 
the following constant emerges: 

pa'^ = p,al\ (45) 

the asterisk meaning a reference time, so far unspecified. 
If A = we have the solution 



whereas if A > 0, 
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with t\ = 2/(37V^). The initial condition is a{t = 0) = 0. 

Equations ()4fi|l and ()47|1 are in agreement with Binetruy et al. [l^. The solution for the 
case A < was however not given in so we give it here: 



4 



sin h 



K 



2t 

: sm — 



(48) 
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where now tx = 2 / {3-y ^J\\\) . 

It is of interest to compare with four- dimensional cosmology. The discovery that the cos- 
mic expansion is accelerating has replaced the Einstein - de Sitter model with the Friedmann 
- Lemaitre model as the standard four- dimensional model This is a fiat universe with 
cosmological constant A4 and pressure-free matter (7 = 1). According to this model, 



a = p^CL 



3 "4 sinh^ ^ 



(49) 



A4 

with = 2/ v^3A4. Equation ()49|1 is seen to correspond to the first term in Eq. ()47j) in view 
of the relationships 

tA, (50) 
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the two first of which follow from a comparison of Eqs. ()28|) and (jHUj) . 



B. Low energy limit 

If p ^ 0", i.e., if the matter density in the universe is much smaller than the intrinsic 
cosmological constant, the Friedmann equation (f28|) reduces to 



H^ = \ + 3^«:Vp, (51) 
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which is is seen to reduce to the four-dimensional Friedmann equation, Eq. ()3U|) . in view of 
the relationships given in Eq. (jMH) . In this case the solutions of five-dimensional cosmology 
can thus be found from the conventional four- dimensional theory. 

Taking the time derivative of Eq. (j51|) and using the energy conservation equation (j43|) . 
we get the following equation for the scalar expansion {9 = 3H): 

e{t) + ^-^e'{t)-^K'ace{t)-^-^x = o. (52) 

If A = 0, in view of /t^cr — > 6^1 this equation is the same as Eq. (31) in Ref. Q]- Now assume 
that A > 0. Integration of Eq. yields 

2 , 29-p + y/K 



t = — =ln 



2^-/?-yZ' 

where /? = (7^/(27) and K = 36 A + Inverting this expression we get 



(53) 



= Ip + iv^coth (^7V^t) • (54) 

Near t = the viscosity is not expected to be important. Ignoring (3 we obtain 6 2/{'-)t) 
when t ^ 0. Thus when p <^ a, five-dimensional cosmology predicts the scalar expansion 
to depend on the form of the equation of state near t = 0. This behaviour of 6 is the same 
as found in four- dimensional cosmology. In the latter case, for any value of the curvature 
parameter k, one has a oc t^/^^T) for t — > 0, resulting in if = 2/(37^), equivalent to the 
expression for 6 given above. 



C. Matter-dominated universe, a = 

In conventional four- dimensional cosmology one assumes, as already mentioned, that 
one epoch was dominated by a non-accelerating expansion in which cold matter dominated 
the universe. This is the Einstein - de Sitter model: a fiat universe with no cosmological 
constant. 

In five- dimensional cosmology the Friedmann equation (j42p in a matter-dominated uni- 
verse would describe a universe in which the term dominates over the ap term. We will in 
the present section simply assume that the brane tension a is zero. Note that the effective 
cosmological constant A of Eq. ()29|) can still dominate over the p^ term if A is large enough, 
thus possibly describing an inflationary universe. 
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The Einstein - de Sitter model corresponded to a vanishing pressure p, i.e., 7 = 1. Here, 
we will admit that there exists a matter pressure, so that 7 will be kept unspecified. As we 
are still outside the radiation epoch we assume, as before, that C = 0. 

The Friedmann equation becomes now 

1 



6b 



(55) 



which, together with the energy conservation equation (|43|) . yields the following equation 
for the scalar expansion: 
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(56) 



We first consider the de Sitter case, A > 0, which implies that also A > in view of the 
assumption a = 0. With the substitution 6 = SV^coshx we obtain 

dx 



(57) 



sinh X — 6 cosh x 

with 6 = K^C/(27). Integrating this equation we obtain a relation between t and x. The 
resulting expression for 6 will depend on the magnitude of the viscosity. Assume first that 
< 27. We then obtain 

'3 



3 
2 



= -v^jitanh f ^-fVXpt 



ficoth ( ^-fVxpt 



where 



B 



/27 + /€2C 



N-(f; 



/ 27 - K^C 

The form of Eq. ()58|) implies that 9 B/{'yf3t) when t — > 0. 
If > 27 we get 



(58) 



(59) 



= ^ { i tan Q7V^ /3 1) + 5 cot (^7V^^ t) } , 
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For an AdS universe, A < implying A < 0, we obtain by an analogous calculation 
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tan -7V^/9t + 5cot -'jV^pt 



(60) 



(61) 



(62) 
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when < 27, and 

3 ^ r 1 /3 _ \ ^ /a 



e = -v^|--^tanh i^-^V^pt^ + 5coth (^-7v^/5t)| , (63) 

when > 27. To our knowledge, these expressions have not been given before. 

Once 6{t) is known, the energy density p{t) is easily obtained from Eq. ()55|) . The scale 
factor a{t) itself is obtained by integrating the equations above. Let us give a{t) for the case 
of a dS space with < 27: 

In the nonviscous case corresponding to 5 = /3 = 1, this expression reduces to 



. , 2t 2 



a 



3-7 — r, ^^'y 



= (65) 

in agreement with Eq. ()47|) (in which the first term is zero). 

It is natural to ask: Which of the above options are most likely to describe real physics in 
the early universe? Most likely one would expect that the bulk viscosity is low, so that the 
inequality < 27 is obeyed. This should support the validity of Eqs. or (jHH). However 
the situation may be more complicated, at least near the Planck time; as we will see more 
closely in subsection D it may happen that the contribution from the bulk viscosity becomes 
quite significant. Moreover, it follows of course mathematically that if 7 happens to be 
close to zero in some epoch in the universe's history, thus corresponding to an approximate 
"vacuum fluid", p = —p, then can easily exceed 27. 

D. Radiation-dominated universe 

In the epoch of the universe's history where radiation dominates, it is no longer permissi- 
ble to neglect the term C/a^ in Friedmann's equation. The equation of state corresponding 
to 7 = 4/3 is 

Ip, (66) 



and the governing equations are 



1/1 1 ,1 n C 



= X+ —K^ap + —K^p^ + ^, (67) 
18 36 a^' ^ ^ 
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P + \pO-ce^ = Q. (68) 

As our primary aim is to study the influence from viscosity, we sliall in tliis subsection ignore 
brane tension and set cr = 0. 

To parametrize tfie viscosity, we make tlie following ansatz for the "kinematic" viscosity 
coefficient C/p associated with C,: 

where a is a constant (the factor 3 is for convenience). This choice is motivated partly by 
analytic tractability of the formalism. In addition, as a more physical motivation, we note 
that it has some similarity with the ansatz made by Murphy three decades ago in his classic 
paper on viscous cosmology 2^. Murphy assumed that C/p = constant. Recalling from 
conventional cosmology that 9 = 3/(2t) when t — > 0, we see that our ansatz describes a 
weaker influence from viscosity near the big bang than does the ansatz of Murphy: Eq. (j69| 
predicts the kinematic viscosity to start from zero at t = and thereafter increase linearly 
with t. 

Inserting Eq. into Eq. we obtain 

pa'-" = P^at"", (70) 

where the asterisk as before means an unspecified instant. The presence of the coefficient 
a > thus implies that p decreases more slowly with increasing a than what is the case for 
a nonviscous fluid. 

The relationship between t and a can be written as an integral {x = a^): 

v^t = 1 (71) 



where 



Qa=[jP.at~-\ . (72) 



We have expressed the integrand in Eq. ()71|) in terms of nondimensional quantities (recall 
that [a] = [a;] = 1, so that [C] = [A] = [Q^] = cm.-^). 

The integral in Eq. (ffT| is easily calculable in three special cases: a = 0,2,4. First, let 
us set a = 0, taking for definiteness A to be positive: 

1 a' + (C/2A) + ^a^ + {C/X)a^ + Qo/X 

vAt=-ln j= . (73) 

4 (C/2A) + ^fo^/X 
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As usual, we have assumed a = for t = 0. For a = 2, we can simply use the same 
expression (fTSj). with the substitutions 

+ Qo^O. (74) 

The effect from viscosity is thus contained in Q2. If a = 4, an analogous substitution can 
be done. This case is however pathological, since it corresponds to p= constant according 
to Eq. (f7n|). 

If A = 0, we have 

t=ir' . (75) 

giving for a = 

t = ^ ijCa' + Qo- V^) . (76) 

The viscous result for a = 2 is obtained by means of the same substitution as above, Eq. ()74p . 

For general values of a the integral in Eq. (f7T|) is of course easily done numerically. Most 
likely, the viscous effects that we consider here are to be attributed to the very early times 
in the history of the universe. As a rather extreme physical example, let us consider the 
evolution of the universe from t = up to the Planck time, Tp = Lp = \fGl = 5.4 x 10~^^ s 
(or 1.6 X 10~^^ cm), and let us identify the asterisk quantities in Eq. (fTOj) with this particular 
instant. Normalizing the nondimensional scale factor such that a^, = 1, meaning that the 
value of a{t) is given in Planck lengths, we may estimate the magnitude of the left hand 
side of Eq. (ffT|l (with t = t^. = Tp) as follows: The magnitude of A, equal to A/6 in our case 



since cr = 0, is taken to be of the same order as the RS fine-tuning energy scale k = y |A|/6 
(ignoring here the minus sign for A). Thus -\/A ~ k. Moreover, we adopt the estimate 
M5 ~ M4 ~ k which is customary in the RS context j21 1 . As M4 = G4 ^^"^ we then get 
\fXTp ~ 1. Moreover, Qa = Qo in view of a* = 1. As very little seems to be known about 
the magnitudes of C/\ and Qo/A we shall here simply assume, as a working hypothesis, 
that these quantities are equal to unity. Thus we obtain from Eq. (f7T|) . when omitting the 
pref actor 1/4, 

1 ~ , "^"^ , (77) 
-'0 Vx^ + x + 

as a condition for determining the magnitude of a. Due to the crudeness of the argument 
we cannot determine a quantitatively from this condition, but we find it of interest to note 
that Eq. (|77j) can be satisfied with a being of the order of unity. 

15 



Finally, let us estimate the magnitude of the viscosity coefficient. We have = SttGs = 
87r/M| ~ Sn/M'i ~ 10"^^ cm^ which, together with the estimate Qo/\ = (^p^)Va ~ 1, 
leads to ~ 10^^° cm^^ (or 10^^ g/cm^). With 6 = 3/(2t) we then obtain at Planck time, 
assuming a ~ 1, 

— ~ 10"^^ cm, C* ~ 10^^ cm-^ (78) 
p* 

Thus the bulk kinematic viscosity C*/p* is small, and the bulk viscosity C* itself is large. 
Moreover, we get k^^^, ~ 1. Somewhat remarkably, this indicates that we are in the transition 
region between the two C < 27 and > 27 considered in subsection C, for the 

matter-dominated universe (cr = 0). There seems to be no simple reason why this should be 
so. 



V. CONCLUDING REMARKS 



We have studied one single flat brane situated at ?/ = in a five-dimensional Randall- 
Sundrum setting, requiring the spatial curvature parameter k to be zero. The isotropic 
fiuid assumed to be present on the brane was assumed to satisfy the conventional equation 
of state p = (7 — l)p, with 7 a constant. The essential new element in our analysis was 
the allowance of a constant bulk viscosity ( in the fiuid. The five-dimensional bulk itself, 
endowed with a five- dimensional cosmological constant A, was assumed not to contain any 
fiuid (or other field). 

The governing equations are the Friedmann equation (f^ . in which A can be looked 
upon as the effective four-dimensional cosmological constant in the five- dimensional theory, 
together with the energy conservation equation (jlUI) . 

In Sec. IV, the main section of our paper, these governing equations were solved under 
various approximations. The nonviscous theory, briefly considered in subsection A, repro- 
duced essentially results obtained earlier SlNl. The results in subsection B in the limit of 



low energy, are formally the same as those found earlier in a four- dimensional context [7(. 
In subsection C, the expressions given analytically for the scalar expansion 6, are to our 
knowledge new results. In subsection D we introduced the ansatz of Eq. (1691) for the bulk 
kinematic viscosity (being related to the ansatz made earlier by Murphy j20|), and derived 
on the basis of this the integral expression in Eq. 1)7111 relating cosmological time t to the 
scale factor a. Some simple analytical results were calculated. Moreover, going to the ex- 
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treme Planck times, our numerical estimates indicated that the viscosity concept may play 
a physical role in the very early universe. 

We emphasize that our theory is built upon the assumption of no flux of energy in the 
fifth direction, as expressed by Eq. (j22j) . 
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